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Abstract. We investigate the coherence of the twin kilohertz quasi-periodic oscillations (kHz QPOs) in the low-mass X-ray 
binary (LMXB) theoretically. The profile of upper kHz QPO, interpreted as Keplerian frequency, is ascribed to the radial extent 
of the kHz QPO emission region, associated with the transitional layer at the magnetosphere-disk boundary, which corresponds 
to the coherence of upper kHz QPO. The theoretical model for Q-factor of upper kHz QPO is applied to the observational data 
of five Atoll and five Z sources, and the consistence is implied. 
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1. Introduction 



, The launch of the Rossi X-ray Timing Explorer (RXTE) led 
to the discovery of Kilohertz quasi-periodic oscillations (kHz 

, QPOs) of X-ray spectra in low-mass X-ray binaries (LMXBs), 
i.e., narrow features in their power density spectra (PDS) (van 
der Klis 2000, 2006). These frequencies, in the range of 200 ~ 
1300 Hz, are the same order of the dynamical time-scales of 
the innermost regions of the accretion flow around the stel- 
lar mass compact objects (van der Klis 2006, 2008), which 
may carry imprints of strong field general relativity phenom- 

' ena (e.g. Kluzniak, Michelson, Wagoner 1990; Miller, Lamb 

; & Psaltis 1998; Kluzniak 1998; Abramowicz et al. 2003). 

' Owing to the expected links with the orbital motion, most 
works about the kHz QPOs focus on the explanation for the 
emission of these frequencies of the orbital Keplerian mo- 
tion (e.g. Miller, Lamb & Psaltis 1998; Stella & Vietri 1998, 
1999; Kluzniak & Abramowicz 2001; Abramowicz et al. 2003; 
Zhang 2004). In addition, the kHz QPOs were found to occur 
in twin peaks usually (upper V2 and lower vi frequency). They 
behave in a rather regular way and follow the tight correlations 
between their frequencies and other observed characteristic fre- 
quencies (see, e.g. Psaltis et al. 1998, 1999a; Psaltis, Belloni & 
van der Klis 1999b; Stella, Vietri & Morsink 1999; Belloni, 
Psaltis & van der Klis 2002; Titarchuk & Wood 2002; Mendez 
& van der KHs 1999, 2000; Mendez et al. 2001; Yu W. F., van 
der Klis M., Jonker P. G. 2001; Yu W. E, van der Klis M. 
2002). Moreover, the correlation between the upper frequency 
and lower frequency across different sources can be roughly fit- 
ted by a power-law function (see e.g. Psaltis et al. 1998, 1999a; 



Zhang et al. 2006a), and also by a linear model (see Belloni, 
Mendez & Homan 2005, 2007). 

The kHz QPOs in LMXB are typical time variable sig- 
nals and peaks with some width in PDS, with the partic- 
ular case of the sharp coherent pulsation 401 Hz and a 
near 401 Hz X-ray burst oscillation frequency are found in 
SAX J1808.4-3658 (Chakrabarty et al. 2003; Wijnands et al. 
2003; Wijnands 2005), whose profiles can be described by a 
Lorentzian Function (see the lower panel in Fig.lTJ, 

Py^Aow/[(y-vof + (w/2f], (1) 

where vo is the peak frequency, w is the full width at half- 
maximum (FWHM), and Aq is the amplitude of this signal. The 
ratio of these two quantities are the quality factor, i.e., 

Q - ^. (2) 

w 

This function is characterized by three characteristic quantities 
which are used to describe the profile of these signals, i.e. cen- 
troid frequency (i.e. peak frequency vq), quality factor (Q = vo / 
FWHM) and the fractional root-mean-squared (rms). The qual- 
ity factor characterizes the coherence of a QPO signal, while 
the rms represents a measure of the signal strength, which is 
proportional to the square root of the peak power contribution 
to the PDS. The large RXTE archive makes possible searches 
for systematic correlations between these three characteristic 
quantities. Using data from RXTE, Barret et al. (2005a) stud- 
ied 4U 1608-52 and revealed a positive correlation between 
lower frequency and its quality factor, up to a maximum of 
about Q ~ 200. Motivated by this idea, BaiTet, Olive & Miller 
(2005b, 2006) studied, in a systematic way, the QPO proper- 
ties of 4U 1636-536 and the dependency of quality factor and 
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Fig. 1. The illustration of Lorentzian function as described in 
Eq. (1) is associated with the sketch map of the magnetosphere- 
disk transition layer. The labels Aq, vq and w are noted in 
Lorentzian function Eq.(l). In the upper panel of figure, the 
stellar radius R, magnetosphere radius r and magnetosphere- 
disk transition layer width 6r are noted, and the kHz QPO pro- 
file and its corresponding positions in the transition layer are 
presented. 



rms on frequency. It is shown that quality factors for the lower 
and upper kHz QPOs of 4U 1636-536 follow different tracks 
in a Q versus frequency plot, i.e. quality factor for the lower 
kHz QPO increases with frequency up to 850 Hz (Q ~ 200) 
and drops precipitously to the highest detected frequencies ~ 
920 Hz (Q ~ 50), while that of the upper kHz QPO increases 
steadily all the way to the highest detectable QPO frequency. 
Moreover, quality factor of the lower QPO is higher than that 
of the upper (Barret, Olive & Miller 2005b,c; 2006). In addi- 
tion, the fractional rms amplitudes of both the upper and lower 
kHz QPOs increase then decrease steadily towards higher fre- 
quencies, with a ceiling in lower kHz QPO (BaiTet, Olive & 
Miller 2005b,c; 2006). A similar behavior was seen from 4U 
1608-52 (Barret et al. 2005a). The rough similarity also was 
extended to 4U 1735-44, 4U 1728-34 (Barret, Olive & Miller 
2006; Boutelier, Barret & Miller 2009; Mendez 2006; Torok 
2009). The complete list of references on QPOs from these 
sources is available in van der Klis (2006). Then Torok (2009) 
studied the difference in rms amplitude between the upper and 
lower kHz QPOs as a function of the frequency ratio (v2/vi). 
They found that the rms amplitudes of the twin peaks become 
equal when the frequencies of the oscillations pass through a 
certain ratio (V2/V1), which is roughly the same for each of 



the sources. It is also predicted that in a more general context, 
the behaviour of the amplitude difference suggests a possible 
energy interchange between the upper and lower QPO modes 
(Torok 2009). 

Theoretically, the quality factor and rms amplitude have 
also been studied but not yet satisfactorily settled. In the past 
years, several works, focussing on the quality factor and rms 
as separated functions of frequencies, have been discussed on 
these properties, and possible consequences for various QPO 
models have been outlined (see, e.g.. Barret et al. 2005a; 
Mendez 2006; Barret, OUve & Miller 2006, for further infor- 
mation and references). Mendez (2006) considered the rela- 
tion between the innermost stable circular orbit (ISCO) and 
the drop of QPO coherence and rms. BaiTet, Olive & Miller 
(2006) discussed the implications of their results and showed 
how the high-frequency drop-off in quality factor can be ac- 
commodated quantitatively in a toy model based on the ap- 
proach to the ISCO. In their toy model, the change of qual- 
ity factor was ascribed to three basic parts: (1) the finite extent 
in radial direction(Aro,i,), (2) the radial drift (Ardrift) during the 
lifetime of the oscillation, (3) the finite time itself (Barret, Olive 
& Miller 2006). It was also proposed that a drop in the ampli- 
tude and quality factor of the QPOs at some limiting frequency 
is a possible signature of the ISCO (Miller, Lamb & Psaltis 
1998; Barret, Olive & Miller 2005b,c). 

The abrupt rise of X-ray flux, as a cause of QPO, may be 
related to the abrupt transformation of environment in accre- 
tion process, which is ascribed to the interaction of magne- 
tosphere dominated regime and gravitational induced energy- 
momentum dominated disk. In this regards, the interpreta- 
tion of the twin kHz QPOs has been ascribed to the orbital 
Keplerian frequency at magnetosphere-disk boundary for up- 
per kHz QPO and MHD Alfven wave propagation frequency 
there for lower kHz QPO, respectively (Zhang 2004; Zhang et 
al. 2007). The quality factor can be attributed to the radial ex- 
tent of the transition zone where kHz QPOs produce (see Fig. 
[T]). In this letter, we investigate how this radial extent leads to 
the profile for upper kHz QPO, which is described in Section 2, 
where the comparison between theory and observations is also 
presented. Section 3 contains the conclusions and discussions. 

2. Formation of High Q Factor and Its Variation 
with Frequency 

According to the MHD Alfven wave oscillation model of kHz 
QPO, the MHD turbulence by the shear flow in the accretion 
disk (e.g., Ruediger & Pipin 2000) will trigger the strong varia- 
tion of plasma energy density and ignite the shear Alfven wave 
motion along the orbit with the loop length of circumference 
27rr at the Keplerian disk orbit radius r (Zhang et al. 2007). 
In the process of NS accretion, there is a transition from the 
spherical accretion with high mass density to a polar cap ac- 
cretion with low mass density (Zhang 2004). It is assumed that 
this transition occuiTed at a certain radius which is called pre- 
feiTed radius, where a MHD tube loop may be formed to con- 
duct the accreted matter to the polar cap of star This critical 
transition may give rise to MHD turbulence and perturb the 
field lines, thus excite Alfven wave oscillation. In such a see- 
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nario, the lower and upper frequency correspond to the AWOF 
with the spherical accretion mass density which coincides with 
the Keplerian orbital frequency and the AWOF with the po- 
lar accretion mass density, respectively. At the preferred ra- 
dius r which is defined by the magnetic pressure matching the 
ram pressure, these two frequencies read (Zhang 2004), respec- 
tively. 



V2 



GM 



1850(i/z)AX5, 



X. 



(3) 
(4) 



Here, the NS mass M is in units of solar mass, and va V^" 
is the AWOF where the area S representing the spherical area 
Sr - AnR^ or the polar cap area Sp - 47rR^(l - Vl - X), re- 
spectively. X = ^ is the ratio between star radius R and disk 
radius r. A = (m/R^^)'^^ with Re = R/10^(cm) and m the mass 
M in the units of solar masses. 

We ascribe the upper kHz QPO frequency profile to the X- 
ray flux strong variations from transition layer, as shown in the 
upper panel of Fig.©. Thus, the range of the upper frequency 
i6v) should be the indication of transition layer radial extent (6r 
or in terms of the scaled quantity X). 



36X l,5r 

OVo = V2 = V2 . 

2X 2r 



(5) 



If we set the QPO FWHM w as corresponding to the scale of 
transition layer width, i.e. w ~ 6v, then the quality factor for 
upper frequency (hereafter upper Q-factor) can be written as. 



22 = ^ = --. 

6v2 3 6r 



(6) 



We consider a transition layer (see upper panel of Fig. ([T])) 
between the innermost Keplerian orbit and the magnetosphere 
(Eisner & Lamb 1977; Naso & Miller 2010ab), in which the 
transition for radial velocity of accretion flow from a Keplerian 
to a corotation with the NS may occur (Titarchuk, Lapidus 
& Mushmov 1999 hereafter TLM; Titarchuk & Osherovich 
1999). But in this region the formation of kinks and shocks due 
to the inhomogeneous density of flow and supersonic motions 
of accretion flow may be responsible for a super-Keplerian ro- 
tation (TLM). After a self-adjustment which is the function 
of this layer, the coupling of the sub-Keplerian flow with the 
super-Keplerian rotation still leads to a Keplerian zone which 
contributes to the radial drift. When the accreted matter fall into 
the innermost Keplerian orbit and transition zone, these plasma 
may strike the magnetospheric boundary and bend the field 
lines, which lead to the change of the magnetospheric shape 
due to the strong ram pressure and can excited the Alfven wave 
oscillation. In the meantime, the plasma is threaded by the field 
lines and fall onto the polar cap. However, on the magneto- 
spheric boundary, the accreted matter carries different veloci- 
ties and has different density from the former arrivals, which 
bring into some instabilities, such as Kelvin-Helmholtz in- 
stability, Rayleigh-Taylor instability, and the inhomogeneities 
caused by the azimuthal component of the field trapped inside 
the inner regions of the disk (Romanova, Kulkarni & Lovelace 



2007). The Kelvin-Helmholtz instability seems to be less sig- 
nificant (Rastatter & Schindler 1999), and we just focus on 
the Rayleigh-Taylor instability (Kulkarni & Romanova 2008). 
The unstable disc-magnetosphere boundary results in the pen- 
etration of the magnetosphere by the disc matter in kinds of 
forms. In addition, the disc matter treading the field lines into 
the magnetosphere from polar cap may move to the equatorial 
region and swell the magnetosphere. On the grounds of this 
idea and the geometry of accretion disc boundary layer (Regev 
& Hougerat 1988), the radial extent direction can be written as 
(Eisner & Lamb 1977; Naso & Miller 2010ab), 



6r — a(r — R), 



(7) 



where the quantity (r-R) stands for the scale of the 
magnetosphere-disk boundary to stelar surface, and a < 1 is 
a constant ratio coefficient. Therefore, we can obtain. 



Q 



-d 



R^ 
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(8) 



3a r 3a 

Using the definition of X and Eq. (O, we can write Q2 with 
respect to the upper frequency, 

V2 



(9) 



1850A' 

For the detected twin kHz QPOs, the mass density param- 
eter A is found to be about 0.7 (e.g. Sco X-1) (Zhang 2004; 
Zhang et al. 2007), except for the two unusual X-ray mil- 
Hsecond pulsar cases SAX J 1804.5-3654 and XTE J 1807-294, 
A=0.45 and 0.4 respectively (Zhang et al. 2010). In most cases 
(except Cir X- 1 ), the position parameter X=R/r lies in the range 
from 0.7 to 0.92, or the kHz QPO emission position radius is 
fromr=l.lR to r=1.4R (Zhang et al. 2010). Averagely, we have 
the following expression to evaluate the upper Q-factor, 

V2 



Q 



(10) 



1300(Hz)' 

To make a comparison between the model and the detec- 
tions of the upper Q-factor, we firstly put the sources whose 
spectrum of kHz QPOs displayed the high coherence together 
in Fig.m We choose the parameter a - 0.35,0.55,0.8, re- 
spectively, and plot the theoretical curves. Then we fit the 
theoretical expression for Q2 (Eq. dlOt ) to the observational 
data and find that the model curve with a - 0.551 fits the 
most data well, i.e. the radial extent of the preferred radius 
is about half of the thickness of the magnetosphere-disk 
boundary to stellar surface (see Fig. |2). It is can be seen 
that five Z sources and five Atoll sources (Hasinger 1989, 
1990) show different coherence trends, which are described 
in Fig.|3l as the same meanings of Fig. |2] but with the sep- 
arate sources for a clearly presentation. The fitting results 
are fisted in Table 1. For the Atoll source 4U 1728-34, there is 
an abrupt drop at v ~ 1050 Hz in the Q-v plot, then it is argued 
as an effect of innermost stable orbit (Barret, Olive & Miller 
2005b, c; 2006) to arise this abrupt kHz QPO profile changes. 

3. Conclusions and Discussions 

We try to ascribe the profile of upper kHz QPO frequency to 
the radial extent of its emission region at magnetosphere-disk 
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Fig. 2. Plot for the upper Q-factor versus its frequency. The data 
are provided by D. Barret and M. Mendez, which have been ex- 
ploited and discussed in the references (BaiTet et al. 2005abc, 
2006, 2007, 2008; Boutelier et al. 2009, 2010; Mendez 2006). 
The theoretical curves of upper Q-factor Q2 of Eq. ( fTOb are 
plotted with the parameter a - 0.35,0.55,0.8 and the aver- 
aged mass density parameter A=0.7 (Zhang 2004; Zhang et al. 
2007). 

Table 1. The fitting results for the observation data. 



Source 


a 


error(Q-) 


X^IDoF 


&uoy{x^ 1 DoF) 


Total 


0.551 


0.030 


3.729 


0.601 


4U 1728-34 


0.400 


0.033 


2.007 


0.452 


4U 1735-44 


0.400 


0.027 


14.095 


5.337 


Sco X-1 


0.650 


0.007 


1.391 


0.539 


Cyg X-2 


0.530 


0.039 


2.186 


1.378 


GX 174-2 


0.800 


0.0340 


6.015 


1.463 


4U 1820-30 


0.530 


0.040 


1.741 


0.522 


GX 340-1-0 


0.600 


0.042 


6.416 


2.068 


GX 5-1 


0.800 


0.022 


1.669 


0.047 


4U 1636-53 


0.530 


0.027 


4.540 


0.450 


4U 1608-52 


0.650 


0.027 


5.183 


1.020 



transition layer, and investigate the evolutions of its quality fac- 
tor as a function of upper frequency. The central frequency 
of upper kHz QPO emits from a preferred radius r, at which 
the magnetic pressure matches the ram pressure of the inward 
disk material, and the radial extent of the prefeiTed radius con- 
tributes to the frequency extension of the upper kHz QPOs. The 
naiTower this extent, the higher coherence of QPO (or upper Q- 
factor). After making the comparisons between the model and 
detected Q-factor data (Fig.|2]i, we find that the radial extent 6r 
accounts for about half of the transition zone, extending from 
the stellar surface to the magnetosphere-disk boundary, e.g. (r- 
R). This means that X-ray flux from this radial extent primar- 
ily contributes to the QPO peak profile (see illustration figure 
FigID- On the quantitatively discussion of the kHz QPO emis- 
sion region, we refer to the recent result by Zhang (2010) that 
the kHz QPOs of most sources emit from the positions at r ~ 20 
km if stellar radius R= 1 5 km is assumed (excluding Cir X- 1 , its 
emission position is slightly far from star, e.g. r ~ 30 km), thus 
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Fig. 3. The same meaning with Fig.|2]but for separated sources 
(five Atoll and five Z sources). The model curves are plotted 
with the different extent parameter a and averaged mass density 
parameter A that were derived from the detecting data (Zhang 
et al. 2007). 



the radial extent of contributing to upper frequency is about 3 
km. 

While making comparisons of model to the Q-factors of 
five Atoll and five Z sources separately, as shown in Fig. [3] 
we notice that there exist a big discrepancy in 4U 1728-34 at 
high frequency, more errors in Cyg X-2 and GX 340H-0. On 
the unusual abrupt drop of upper Q-factor of 4U 1728-34, it is 
ascribed as the effect of ISCO (BaiTet, Olive & Miller 2006). 
On the upper Q-factors of Cyg X-2 and GX 340 + 0, we argue 
that there might not only exists a radial extent but also a verti- 
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cal extent to contribute to the QPO emission regions, and if the 
thickness extent of transition layer is considered it might arise 
new ingredients on the QPO profiles. On the thickness of accre- 
tion disk and QPO phenomenon, we refer to the recent work 
by Chakrabarti et al. (2009). We think that the complete de- 
scription of QPO profile should take both radial and thickness 
extents around preferred radius into account, since the pileup 
blob of accreting plasma interacts with magnetosphere to pro- 
duce the fruitful X-rays in a spacial domain. Moreover, it is 
remarked that our model on the profile of kHz QPO is concen- 
trated to the upper frequency Q2, which cannot be automati- 
cally appUed to the lower frequency Q\, and in fact the obser- 
vations of Q2 and Qi are apparently different. We ascribe the 
profile width of V2 to the radial width of the transition layer of 
magnetosphere-disk, which is scaled by 5r - r-R, but the pro- 
file width of v\ should be not related to this width. Therefore, 
there does not exist a relation between Qi and Q2, although a 
relation between the upper and lower frequencies can be clearly 
expressed. It is a subsequent work to figure out the mechanism 
for the profile of lower kHz QPO gi ■ 
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